
R : Copyright 2003, The R Development Core Team
Version 1.8.0  (2003-10-08)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

> invisible(options(echo = TRUE))
> # estimate John Jackson's model using multinomRob
> #
> 
> options(width=130);
> 
> library(multinomRob)
Loading required package: rgenoud 
Loading required package: MASS 
Loading required package: mvtnorm 
> 
> # read John Jackson's replication dataset (PA 2002, Winter)
> #
> jjdat <- read.table(file="pareplic.asc", header=T);
> 
> #> names(jjdat)
> # [1] "Voivodship"    "ID"            "UD.KLD"        "SLD"
> # [5] "PSL"           "UP"            "Catholic"      "Other"
> # [9] "Constant"      "New.Jobs"      "Unemp"         "SOE"
> #[13] "Church.Att"    "Schooling"     "Age"           "Farmer"
> #[17] "Village"       "Suchoka.UD"    "Pawlak.PSL"    "Kwasn.SLD"
> #[21] "Zioqkowska.UP" "Voters"
> 
> jjuse <- jjdat;
> #total number of voters
> jjuse$TotalVotes   <- jjdat$Voters * 1000;
> jjuse$SLD      <- jjdat$SLD/100      *jjuse$TotalVotes;
> jjuse$PSL      <- jjdat$PSL/100      *jjuse$TotalVotes;
> jjuse$UP       <- jjdat$UP/100       *jjuse$TotalVotes;
> jjuse$Catholic <- jjdat$Catholic/100 *jjuse$TotalVotes;
> jjuse$Other    <- jjdat$Other/100    *jjuse$TotalVotes;
> jjuse$UD.KLD   <- jjdat$UD.KLD/100   *jjuse$TotalVotes;
> 
> jjuse$Age  <-  jjdat$Age /10;
> 
> ## Set parameters for Genoud
> zz.genoud.parms <- list( pop.size             = 3000,
+                         wait.generations      = 30,
+                         P5                    = 500,
+                         scale.domains         = 4,
+                         MemoryMatrix          = TRUE
+                         )
> 
> # estimate the model as done in the Mebane-Sekhon AJPS 2004 paper
> 
> load("jj.a.RData");
> starts <- muljj$genoud$par;
> #starts[c(8,18,26,41)] <- muljj$mtanh$coeffvec[c(8,18,26,41)];
> 
> muljj <-
+   multinomRob(
+     list(SLD ~      New.Jobs + Unemp + SOE + Church.Att + Schooling + Age + Kwasn.SLD ,
+          PSL ~      New.Jobs + Unemp + SOE + Church.Att + Schooling + Age + Farmer + Village + Pawlak.PSL ,
+          UP ~       New.Jobs + Unemp + SOE + Church.Att + Schooling + Age + Zioqkowska.UP ,
+          Catholic ~ New.Jobs + Unemp + SOE + Church.Att + Schooling + Age ,
+          Other ~    New.Jobs + Unemp + SOE + Church.Att + Schooling + Age ,
+          UD.KLD ~   0 + Suchoka.UD ),
+     jjuse,
+     starting.values = starts,
+     genoud.parms = zz.genoud.parms,
+     print.level = 3);

Your Model (xvec):
                                                                       SLD PSL UP Catholic Other UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD   1   1  1        1     1      1
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                          1   1  1        1     1      0
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                         1   1  1        1     1      0
SOE/SOE/SOE/SOE/SOE/NA                                                   1   1  1        1     1      0
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA                1   1  1        1     1      0
Schooling/Schooling/Schooling/Schooling/Schooling/NA                     1   1  1        1     1      0
Age/Age/Age/Age/Age/NA                                                   1   1  1        1     1      0
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                  1   1  1        0     0      0
NA/Village/NA/NA/NA/NA                                                   0   1  0        0     0      0
NA/Pawlak.PSL/NA/NA/NA/NA                                                0   1  0        0     0      0
multinomRob(): Starting Values 
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.67449612  1.7759016  1.7839568 -6.9921414  3.0286856
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.50118106 -7.6977606 -3.6393789 -3.2561745 -5.1886918
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.66006377  0.8478765  0.7781216  2.2463576  0.5290754
SOE/SOE/SOE/SOE/SOE/NA                                                  0.53772545  0.5102214  0.7312367  0.1792841  0.8100803
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.18235315 -1.2092785 -1.7802717  2.9983352  0.1435616
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38926536 -0.3465239 -0.2813651 -0.3583765 -0.4011777
Age/Age/Age/Age/Age/NA                                                  0.07286005  0.3793209  0.3484988  1.7512541  0.4776720
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.57521745  1.9328898  0.7405775  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.0919025  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9798359  0.0000000  0.0000000  0.0000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD -0.3312504
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                         0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                  0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA               0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                    0.0000000
Age/Age/Age/Age/Age/NA                                                  0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000
multinomRob(): starting fit = 42.20445 


Fri Mar  5 15:05:57 2004
Domains:
 -1.702349e+01   <=  X1   <=    2.837248e+01 
 -3.250591e+01   <=  X2   <=    1.950354e+01 
 -1.980191e+00   <=  X3   <=    3.300319e+00 
 -1.613176e+00   <=  X4   <=    2.688627e+00 
 -1.091177e+01   <=  X5   <=    6.547059e+00 
 -1.946327e+00   <=  X6   <=    1.167796e+00 
 -2.185801e-01   <=  X7   <=    3.643002e-01 
 -1.725652e+00   <=  X8   <=    2.876087e+00 
 -5.327705e+00   <=  X9   <=    8.879508e+00 
 -3.848880e+01   <=  X10  <=    2.309328e+01 
 -2.543630e+00   <=  X11  <=    4.239383e+00 
 -1.530664e+00   <=  X12  <=    2.551107e+00 
 -6.046392e+00   <=  X13  <=    3.627835e+00 
 -1.732620e+00   <=  X14  <=    1.039572e+00 
 -1.137963e+00   <=  X15  <=    1.896604e+00 
 -5.798670e+00   <=  X16  <=    9.664449e+00 
 -6.275708e+00   <=  X17  <=    1.045951e+01 
 -2.939508e+00   <=  X18  <=    4.899180e+00 
 -5.351870e+00   <=  X19  <=    8.919784e+00 
 -1.819689e+01   <=  X20  <=    1.091814e+01 
 -2.334365e+00   <=  X21  <=    3.890608e+00 
 -2.193710e+00   <=  X22  <=    3.656183e+00 
 -8.901359e+00   <=  X23  <=    5.340815e+00 
 -1.406825e+00   <=  X24  <=    8.440952e-01 
 -1.045496e+00   <=  X25  <=    1.742494e+00 
 -2.221733e+00   <=  X26  <=    3.702888e+00 
 -3.496071e+01   <=  X27  <=    2.097642e+01 
 -1.628087e+01   <=  X28  <=    9.768524e+00 
 -6.739073e+00   <=  X29  <=    1.123179e+01 
 -5.378524e-01   <=  X30  <=    8.964206e-01 
 -8.995006e+00   <=  X31  <=    1.499168e+01 
 -1.791882e+00   <=  X32  <=    1.075129e+00 
 -5.253762e+00   <=  X33  <=    8.756270e+00 
 -9.086057e+00   <=  X34  <=    1.514343e+01 
 -2.594346e+01   <=  X35  <=    1.556608e+01 
 -1.587226e+00   <=  X36  <=    2.645377e+00 
 -2.430241e+00   <=  X37  <=    4.050401e+00 
 -4.306849e-01   <=  X38  <=    7.178081e-01 
 -2.005889e+00   <=  X39  <=    1.203533e+00 
 -1.433016e+00   <=  X40  <=    2.388360e+00 
 -1.656252e+00   <=  X41  <=    9.937511e-01 

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	179
	(2) Uniform Mutation.................. 	176
	(3) Boundary Mutation................. 	176
	(4) Non-Uniform Mutation.............. 	176
	(5) Polytope Crossover................ 	1764
	(6) Multiple Point Simple Crossover... 	176
	(7) Whole Non-Uniform Mutation........ 	176
	(8) Heuristic Crossover............... 	176
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 30
Population size       : 3000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
 5.67	-6.50	 0.66	 0.54	-2.18	-0.39	 0.07	 0.58	 1.78	-7.70	 0.85	 0.51	-1.21	-0.35	 0.38	 1.93	 2.09	 0.98	 1.78	-3.64	 0.78	 0.73	-1.78	-0.28	 0.35	 0.74	-6.99	-3.26	 2.25	 0.18	 3.00	-0.36	 1.75	 3.03	-5.19	 0.53	 0.81	 0.14	-0.40	 0.48	-0.33	
fitness = 4.220445e+01

 2.64	-22.23	-0.73	-0.73	-9.01	 0.28	 0.35	 0.29	-1.56	-36.02	 3.94	 0.78	 2.84	-1.60	 1.55	 8.42	-5.13	 0.36	 3.02	 0.98	-1.70	 1.59	 0.57	-1.15	 1.62	 2.08	-16.01	 6.49	-3.29	-0.13	 9.44	 0.18	 1.08	 6.89	-19.76	 0.18	 2.23	 0.25	-1.30	 1.39	 0.53	
fitness = 5.509317e+02


The worst fit of the population is: 9.592898e+17




Generation#	    Solution Value

      0 	4.220445e+01
XXX
XXX

Soft Generation Wait Limit Hit.
No Improvement in 30 Generations

Best Fit Found at Generation 1
Fit Value = 4.220445e+01


Parameters at the Solution (value, gradient):

 X[ 1] :	5.674496e+00	G[ 1] :	1.381915e+02
 X[ 2] :	-6.501181e+00	G[ 2] :	1.828521e+01
 X[ 3] :	6.600638e-01	G[ 3] :	-2.040970e+01
 X[ 4] :	5.377255e-01	G[ 4] :	4.024517e+01
 X[ 5] :	-2.182353e+00	G[ 5] :	1.089258e+02
 X[ 6] :	-3.892654e-01	G[ 6] :	1.543279e+03
 X[ 7] :	7.286005e-02	G[ 7] :	6.235179e+02
 X[ 8] :	5.752175e-01	G[ 8] :	0.000000e+00
 X[ 9] :	1.775902e+00	G[ 9] :	-5.069385e+01
 X[10] :	-7.697761e+00	G[10] :	-8.280534e+00
 X[11] :	8.478765e-01	G[11] :	5.552585e+00
 X[12] :	5.102214e-01	G[12] :	-1.818337e+01
 X[13] :	-1.209278e+00	G[13] :	-3.148948e+01
 X[14] :	-3.465239e-01	G[14] :	-6.848390e+02
 X[15] :	3.793209e-01	G[15] :	-2.286472e+02
 X[16] :	1.932890e+00	G[16] :	1.043626e+01
 X[17] :	2.091903e+00	G[17] :	-2.009374e+00
 X[18] :	9.798359e-01	G[18] :	0.000000e+00
 X[19] :	1.783957e+00	G[19] :	1.613802e+02
 X[20] :	-3.639379e+00	G[20] :	1.504711e+01
 X[21] :	7.781216e-01	G[21] :	3.148174e+01
 X[22] :	7.312367e-01	G[22] :	2.631200e+01
 X[23] :	-1.780272e+00	G[23] :	1.176155e+02
 X[24] :	-2.813651e-01	G[24] :	1.817990e+03
 X[25] :	3.484988e-01	G[25] :	7.236617e+02
 X[26] :	7.405775e-01	G[26] :	-5.972189e+01
 X[27] :	-6.992141e+00	G[27] :	1.542489e+00
 X[28] :	-3.256175e+00	G[28] :	6.066867e-01
 X[29] :	2.246358e+00	G[29] :	-1.266417e+00
 X[30] :	1.792841e-01	G[30] :	1.329285e+00
 X[31] :	2.998335e+00	G[31] :	9.691880e-01
 X[32] :	-3.583765e-01	G[32] :	3.068130e+01
 X[33] :	1.751254e+00	G[33] :	7.577400e+00
 X[34] :	3.028686e+00	G[34] :	-5.571936e+02
 X[35] :	-5.188692e+00	G[35] :	-6.493638e+01
 X[36] :	5.290754e-01	G[36] :	-3.928630e+01
 X[37] :	8.100803e-01	G[37] :	-3.153916e+02
 X[38] :	1.435616e-01	G[38] :	-2.805744e+02
 X[39] :	-4.011777e-01	G[39] :	1.364850e+02
 X[40] :	4.776720e-01	G[40] :	-1.104469e+01
 X[41] :	-3.318405e-01	G[41] :	-5.513578e+01

Solution Found Generation 1
Number of Generations Run 31


Fri Mar  5 15:08:37 2004
Total run time : 0 hours 2 minutes and 40 seconds


LQD Results:
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.67449612  1.7759016  1.7839568 -6.9921414  3.0286856
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.50118106 -7.6977606 -3.6393789 -3.2561745 -5.1886918
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.66006377  0.8478765  0.7781216  2.2463576  0.5290754
SOE/SOE/SOE/SOE/SOE/NA                                                  0.53772545  0.5102214  0.7312367  0.1792841  0.8100803
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.18235315 -1.2092785 -1.7802717  2.9983352  0.1435616
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38926536 -0.3465239 -0.2813651 -0.3583765 -0.4011777
Age/Age/Age/Age/Age/NA                                                  0.07286005  0.3793209  0.3484988  1.7512541  0.4776720
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.57521745  1.9328898  0.7405775  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.0919025  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9798359  0.0000000  0.0000000  0.0000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD -0.3318405
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                         0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                  0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA               0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                    0.0000000
Age/Age/Age/Age/Age/NA                                                  0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000

LQD sigma: 42.20445 


(multinomTanh):
Tanh Estimates
                                                                              SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.6568082  1.7407704  1.7950218 -6.8295048  3.0699272
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.5345051 -7.7241042 -3.6312504 -3.0662426 -5.1688484
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.6875825  0.9049460  0.7924629  2.3348648  0.4893721
SOE/SOE/SOE/SOE/SOE/NA                                                  0.5468140  0.5247310  0.7350656  0.3097461  0.7967586
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.1491537 -1.2446122 -1.7759156  2.9653124  0.1354074
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.3892047 -0.3462772 -0.2809889 -0.3626082 -0.4014914
Age/Age/Age/Age/Age/NA                                                  0.0720115  0.3871414  0.3445385  1.7205414  0.4721936
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.5694350  1.9411410  0.7348873  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000  2.1061750  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000  0.9728579  0.0000000  0.0000000  0.0000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD -0.3309837
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                         0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                  0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA               0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                    0.0000000
Age/Age/Age/Age/Age/NA                                                  0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000

Tanh Sandwich SEs
                                                                              SLD        PSL         UP  Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 2.04718715 3.56107339 2.33423515 5.0766439 2.01471096
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        2.22110914 3.25306895 1.41092916 3.2814963 1.68331922
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.80840946 1.22402880 0.79669677 1.6316265 0.91986465
SOE/SOE/SOE/SOE/SOE/NA                                                 0.32231960 0.38683369 0.24454628 0.8163544 0.30148879
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.39479683 0.50362956 0.42625546 0.6369009 0.40341748
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.08938293 0.11064989 0.07495529 0.1604958 0.08642515
Age/Age/Age/Age/Age/NA                                                 0.34598267 0.61821180 0.35859573 0.8282866 0.31817216
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.05179565 0.93140939 0.07916530 0.0000000 0.00000000
NA/Village/NA/NA/NA/NA                                                 0.00000000 0.81211954 0.00000000 0.0000000 0.00000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.00000000 0.07158772 0.00000000 0.0000000 0.00000000
                                                                          UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 0.2369377
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                 0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.0000000
Age/Age/Age/Age/Age/NA                                                 0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.0000000
NA/Village/NA/NA/NA/NA                                                 0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.0000000

TANH sigma: 35.55004 


**********************************

Iteration: 1 
using best non-genoud LQD sigma. Provided by Tanh.


Fri Mar  5 15:08:39 2004
Domains:
 -5.656808e+00   <=  X1   <=    1.697042e+01 
 -1.960352e+01   <=  X2   <=    6.534505e+00 
 -6.875825e-01   <=  X3   <=    2.062748e+00 
 -5.468140e-01   <=  X4   <=    1.640442e+00 
 -6.447461e+00   <=  X5   <=    2.149154e+00 
 -1.167614e+00   <=  X6   <=    3.892047e-01 
 -7.201150e-02   <=  X7   <=    2.160345e-01 
 -5.694350e-01   <=  X8   <=    1.708305e+00 
 -1.740770e+00   <=  X9   <=    5.222311e+00 
 -2.317231e+01   <=  X10  <=    7.724104e+00 
 -9.049460e-01   <=  X11  <=    2.714838e+00 
 -5.247310e-01   <=  X12  <=    1.574193e+00 
 -3.733837e+00   <=  X13  <=    1.244612e+00 
 -1.038832e+00   <=  X14  <=    3.462772e-01 
 -3.871414e-01   <=  X15  <=    1.161424e+00 
 -1.941141e+00   <=  X16  <=    5.823423e+00 
 -2.106175e+00   <=  X17  <=    6.318525e+00 
 -9.728579e-01   <=  X18  <=    2.918574e+00 
 -1.795022e+00   <=  X19  <=    5.385065e+00 
 -1.089375e+01   <=  X20  <=    3.631250e+00 
 -7.924629e-01   <=  X21  <=    2.377389e+00 
 -7.350656e-01   <=  X22  <=    2.205197e+00 
 -5.327747e+00   <=  X23  <=    1.775916e+00 
 -8.429666e-01   <=  X24  <=    2.809889e-01 
 -3.445385e-01   <=  X25  <=    1.033616e+00 
 -7.348873e-01   <=  X26  <=    2.204662e+00 
 -2.048851e+01   <=  X27  <=    6.829505e+00 
 -9.198728e+00   <=  X28  <=    3.066243e+00 
 -2.334865e+00   <=  X29  <=    7.004594e+00 
 -3.097461e-01   <=  X30  <=    9.292382e-01 
 -2.965312e+00   <=  X31  <=    8.895937e+00 
 -1.087825e+00   <=  X32  <=    3.626082e-01 
 -1.720541e+00   <=  X33  <=    5.161624e+00 
 -3.069927e+00   <=  X34  <=    9.209782e+00 
 -1.550655e+01   <=  X35  <=    5.168848e+00 
 -4.893721e-01   <=  X36  <=    1.468116e+00 
 -7.967586e-01   <=  X37  <=    2.390276e+00 
 -1.354074e-01   <=  X38  <=    4.062221e-01 
 -1.204474e+00   <=  X39  <=    4.014914e-01 
 -4.721936e-01   <=  X40  <=    1.416581e+00 
 -9.929510e-01   <=  X41  <=    3.309837e-01 

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	179
	(2) Uniform Mutation.................. 	176
	(3) Boundary Mutation................. 	176
	(4) Non-Uniform Mutation.............. 	176
	(5) Polytope Crossover................ 	1764
	(6) Multiple Point Simple Crossover... 	176
	(7) Whole Non-Uniform Mutation........ 	176
	(8) Heuristic Crossover............... 	176
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 30
Population size       : 3000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
 5.66	-6.53	 0.69	 0.55	-2.15	-0.39	 0.07	 0.57	 1.74	-7.72	 0.90	 0.52	-1.24	-0.35	 0.39	 1.94	 2.11	 0.97	 1.80	-3.63	 0.79	 0.74	-1.78	-0.28	 0.34	 0.73	-6.83	-3.07	 2.33	 0.31	 2.97	-0.36	 1.72	 3.07	-5.17	 0.49	 0.80	 0.14	-0.40	 0.47	-0.33	
fitness = 4.276149e+01

-1.44	 5.70	 1.26	 0.55	-2.98	 0.13	 0.17	-0.02	-0.53	-16.74	 0.35	 1.00	-2.60	 0.03	-0.19	 4.67	 2.70	 0.21	 1.75	-1.16	 0.01	 1.44	-1.85	-0.58	 0.74	-0.09	-11.00	-4.07	 6.78	 0.70	-0.86	-0.65	 3.48	-2.87	-3.25	 0.09	 0.81	 0.32	-0.68	 1.19	-0.45	
fitness = 2.066459e+02


The worst fit of the population is: 1.659253e+10




Generation#	    Solution Value

      0 	4.276149e+01
XXX
XXX
XXX
      1 	4.244742e+01
      2 	4.240380e+01
XXX
XXX
XXX
XXX
      3 	4.236597e+01
      4 	4.236356e+01
      6 	4.236333e+01
XXX
XXX
XXX
XXX
      7 	4.236148e+01
XXX
XXX
XXX
XXX
XXX
      9 	4.236034e+01
     10 	4.235446e+01
XXX
XXX
XXX
XXX
     11 	4.235183e+01
     13 	4.235160e+01
     14 	4.235145e+01
     15 	4.235118e+01
     16 	4.235108e+01
     17 	4.235107e+01
     18 	4.235107e+01
     19 	4.235107e+01
     20 	4.235107e+01
     21 	4.235107e+01
     22 	4.235107e+01
     23 	4.235107e+01
     24 	4.235107e+01
     25 	4.235107e+01
     26 	4.235107e+01
     27 	4.235107e+01
     28 	4.235107e+01
     29 	4.235107e+01
     30 	4.235107e+01
     32 	4.234672e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     33 	4.234331e+01
     34 	4.234271e+01
     35 	4.234266e+01
     36 	4.234209e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     37 	4.234015e+01
     38 	4.233959e+01
     39 	4.233919e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     40 	4.233918e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     41 	4.233901e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     42 	4.233899e+01
     43 	4.233899e+01
     44 	4.233898e+01
     45 	4.233897e+01
     46 	4.233897e+01
     47 	4.233897e+01
     48 	4.233896e+01
     49 	4.233896e+01
     50 	4.233896e+01
     51 	4.233896e+01
     52 	4.233889e+01
     53 	4.233883e+01
     54 	4.233831e+01
     55 	4.233806e+01
     56 	4.233795e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     57 	4.233792e+01
XXX
XXX
     58 	4.233781e+01
     59 	4.233780e+01
XXX
XXX
XXX
     60 	4.233772e+01
     61 	4.233771e+01
     62 	4.233766e+01
     63 	4.233762e+01
XXX
XXX
XXX
     64 	4.233757e+01
XXX
XXX
XXX
XXX
     65 	4.233752e+01
     66 	4.233742e+01
     67 	4.233736e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     68 	4.233730e+01
     69 	4.233728e+01
     70 	4.233726e+01
     71 	4.233696e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     72 	4.233692e+01
XXX
XXX
XXX
XXX
     73 	4.233691e+01
XXX
XXX
     74 	4.233687e+01
     75 	4.233683e+01
XXX
XXX
     76 	4.233675e+01
XXX
XXX
     77 	4.233674e+01
XXX
XXX
XXX
XXX
XXX
XXX
     78 	4.233671e+01
     79 	4.233665e+01
XXX
XXX
     80 	4.233662e+01
     81 	4.233655e+01
     82 	4.233652e+01
XXX
XXX
     83 	4.233643e+01
     84 	4.233641e+01
XXX
XXX
     85 	4.233600e+01
XXX
XXX
XXX
     86 	4.233578e+01
XXX
XXX
     87 	4.233575e+01
     88 	4.233560e+01
     89 	4.233543e+01
     90 	4.233525e+01
XXX
XXX
     91 	4.233516e+01
     92 	4.233512e+01
XXX
XXX
XXX
XXX
     93 	4.233503e+01
     94 	4.233500e+01
     95 	4.233490e+01
XXX
XXX
XXX
     96 	4.233438e+01
XXX
XXX
XXX
XXX
XXX
     97 	4.233379e+01
XXX
XXX
XXX
     98 	4.233364e+01
     99 	4.233287e+01
XXX
XXX
    100 	4.233271e+01

Increasing Soft Maximum Generation Limit by 100 (MaxGenerations) to 200.
I'm doing this because the fitness is still improving.

XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    101 	4.233269e+01
XXX
XXX
XXX
    102 	4.233264e+01
    103 	4.233207e+01
    104 	4.233203e+01
    105 	4.233178e+01
XXX
XXX
XXX
    106 	4.233062e+01
XXX
XXX
XXX
XXX
    107 	4.233045e+01
XXX
XXX
XXX
XXX
    108 	4.232965e+01
    109 	4.232895e+01
    110 	4.232866e+01
    111 	4.232842e+01
    112 	4.232832e+01
    113 	4.232807e+01
    114 	4.232800e+01
    115 	4.232789e+01
    116 	4.232788e+01
    117 	4.232788e+01
    118 	4.232787e+01
    119 	4.232785e+01
    120 	4.232778e+01
    121 	4.232772e+01
    122 	4.232771e+01
    123 	4.232771e+01
    124 	4.232771e+01
    125 	4.232771e+01
    126 	4.232771e+01
    127 	4.232771e+01
    128 	4.232771e+01
    129 	4.232771e+01
    130 	4.232770e+01
    131 	4.232770e+01
    132 	4.232770e+01
    133 	4.232770e+01
    134 	4.232770e+01
    135 	4.232770e+01
    136 	4.232770e+01
    149 	4.232755e+01
    150 	4.232730e+01
    151 	4.232706e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    152 	4.232700e+01
    153 	4.232675e+01
    154 	4.232671e+01
    155 	4.232532e+01
    156 	4.232508e+01
    157 	4.232347e+01
    158 	4.232185e+01
    159 	4.232060e+01
    160 	4.231932e+01
    161 	4.231884e+01
    162 	4.231776e+01
    163 	4.231736e+01
    164 	4.231646e+01

WARNING: Turning Off MemoryMatrix because memory usage is too great.

    165 	4.231494e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    166 	4.231418e+01
    167 	4.231394e+01
    168 	4.231345e+01
    169 	4.231332e+01
    170 	4.231317e+01
    171 	4.231309e+01
    172 	4.231299e+01
    173 	4.231295e+01
    174 	4.231283e+01
    175 	4.231267e+01
    176 	4.231262e+01
    177 	4.231256e+01
    178 	4.231240e+01
    179 	4.231239e+01
    180 	4.231235e+01
    181 	4.231235e+01
    182 	4.231235e+01
    183 	4.231234e+01
    184 	4.231234e+01
    185 	4.231234e+01
    186 	4.231234e+01
    187 	4.231233e+01
    188 	4.231233e+01
    189 	4.231233e+01
    190 	4.231232e+01
    191 	4.231232e+01
    192 	4.231232e+01
    193 	4.231232e+01
    194 	4.231232e+01
    195 	4.231232e+01
    196 	4.231232e+01
    197 	4.231232e+01
    198 	4.231232e+01
    199 	4.231224e+01
XXX
XXX
XXX
XXX
XXX
XXX
    200 	4.231200e+01

Increasing Soft Maximum Generation Limit by 200 (MaxGenerations) to 400.
I'm doing this because the fitness is still improving.

    201 	4.231190e+01
    202 	4.231188e+01
    203 	4.231187e+01
    204 	4.231187e+01
    205 	4.231187e+01
    206 	4.231173e+01
    207 	4.231157e+01
    208 	4.231157e+01
    209 	4.231156e+01
    210 	4.231156e+01
    211 	4.231156e+01
    212 	4.231098e+01
    213 	4.231041e+01
    214 	4.231035e+01
    215 	4.231016e+01
    216 	4.231010e+01
    217 	4.231004e+01
XXX
    218 	4.230967e+01
    219 	4.230964e+01
    220 	4.230963e+01
XXX
XXX
    221 	4.230963e+01
XXX
    222 	4.230963e+01
XXX
    223 	4.230963e+01
    224 	4.230963e+01
    225 	4.230962e+01
    226 	4.230962e+01
    227 	4.230962e+01
    228 	4.230962e+01
    229 	4.230962e+01
    230 	4.230962e+01
    231 	4.230962e+01
    232 	4.230962e+01
    233 	4.230962e+01
    234 	4.230962e+01
    235 	4.230962e+01
    236 	4.230962e+01
    237 	4.230962e+01
    238 	4.230962e+01
    239 	4.230958e+01
    240 	4.230956e+01
    241 	4.230955e+01
    242 	4.230955e+01
XXX
    243 	4.230955e+01
    244 	4.230949e+01
    245 	4.230949e+01
    246 	4.230946e+01
    247 	4.230916e+01
    248 	4.230916e+01
    249 	4.230916e+01
    250 	4.230916e+01
    251 	4.230915e+01
    252 	4.230914e+01
    253 	4.230914e+01
    254 	4.230914e+01
    255 	4.230914e+01
    256 	4.230913e+01
    257 	4.230913e+01
    258 	4.230913e+01
    259 	4.230913e+01
    260 	4.230913e+01
    261 	4.230913e+01
    262 	4.230913e+01
    269 	4.230913e+01
    270 	4.230913e+01
    271 	4.230913e+01
    272 	4.230912e+01
    273 	4.230912e+01
    274 	4.230912e+01
    275 	4.230912e+01
    276 	4.230912e+01
    277 	4.230912e+01
    278 	4.230912e+01
    279 	4.230912e+01
    280 	4.230912e+01
    297 	4.230911e+01
    298 	4.230911e+01
    299 	4.230911e+01
    300 	4.230911e+01
    301 	4.230911e+01
    302 	4.230911e+01
    303 	4.230911e+01
    304 	4.230911e+01
    305 	4.230911e+01
    306 	4.230911e+01
    307 	4.230911e+01
    308 	4.230828e+01
XXX
    309 	4.230828e+01
    310 	4.230819e+01
XXX
    311 	4.230818e+01
    312 	4.230815e+01
    313 	4.230815e+01
    314 	4.230814e+01
    315 	4.230811e+01
    316 	4.230810e+01
    317 	4.230808e+01
    318 	4.230808e+01
    319 	4.230808e+01
    320 	4.230807e+01
    321 	4.230807e+01
    322 	4.230807e+01
    323 	4.230806e+01
    324 	4.230805e+01
    325 	4.230805e+01
    326 	4.230804e+01
    327 	4.230804e+01
    328 	4.230804e+01
    329 	4.230803e+01
    330 	4.230794e+01
    331 	4.230771e+01
    332 	4.230769e+01
    333 	4.230767e+01
    334 	4.230767e+01
    335 	4.230765e+01
    336 	4.230736e+01
    337 	4.230734e+01
    338 	4.230732e+01
    339 	4.230712e+01
    340 	4.230696e+01
XXX
    341 	4.230691e+01
    342 	4.230690e+01
    343 	4.230689e+01
    344 	4.230688e+01
    345 	4.230685e+01
XXX
    346 	4.230684e+01
XXX
    347 	4.230679e+01
    348 	4.230679e+01
XXX
    349 	4.230679e+01
XXX
    350 	4.230678e+01
XXX
    351 	4.230678e+01
XXX
    352 	4.230678e+01
XXX
    353 	4.230678e+01
XXX
    354 	4.230678e+01
XXX
    355 	4.230678e+01
XXX
    356 	4.230678e+01
XXX
    357 	4.230678e+01
XXX
    358 	4.230678e+01
XXX
    359 	4.230678e+01
XXX
    360 	4.230678e+01
XXX
    361 	4.230678e+01
XXX
    362 	4.230678e+01
XXX
    363 	4.230678e+01
XXX
    364 	4.230678e+01
XXX
    365 	4.230678e+01
XXX
    366 	4.230678e+01
XXX
    367 	4.230678e+01
XXX
    368 	4.230678e+01
XXX
    369 	4.230678e+01
XXX
    370 	4.230678e+01
XXX
    371 	4.230678e+01
XXX
    372 	4.230678e+01
XXX
    373 	4.230678e+01
XXX
    374 	4.230678e+01
XXX
    375 	4.230678e+01
XXX
    376 	4.230678e+01
XXX
    377 	4.230678e+01
XXX
    378 	4.230678e+01
XXX
    379 	4.230678e+01
XXX
    380 	4.230678e+01
XXX
    381 	4.230678e+01
XXX
    382 	4.230678e+01
XXX
    383 	4.230678e+01
XXX
    384 	4.230678e+01
XXX
    385 	4.230678e+01
XXX
    386 	4.230678e+01
XXX
    387 	4.230678e+01
XXX
    388 	4.230678e+01
XXX
    389 	4.230678e+01
XXX
    390 	4.230678e+01
XXX
    391 	4.230678e+01
XXX
    392 	4.230678e+01
XXX
    393 	4.230678e+01
XXX
    394 	4.230678e+01
XXX
    395 	4.230678e+01
XXX
    396 	4.230678e+01
XXX
    397 	4.230678e+01
XXX
    398 	4.230678e+01
    399 	4.230678e+01
    400 	4.230678e+01

Increasing Soft Maximum Generation Limit by 400 (MaxGenerations) to 800.
I'm doing this because the fitness is still improving.

    401 	4.230678e+01
    402 	4.230678e+01
    403 	4.230678e+01
    404 	4.230678e+01
    405 	4.230678e+01
    406 	4.230678e+01
    407 	4.230678e+01
    408 	4.230678e+01
    409 	4.230678e+01
    410 	4.230678e+01
    411 	4.230678e+01
    412 	4.230678e+01
    413 	4.230678e+01
    414 	4.230678e+01
    415 	4.230678e+01
    417 	4.230678e+01
    418 	4.230678e+01
    419 	4.230678e+01
    420 	4.230678e+01
    421 	4.230678e+01

Soft Generation Wait Limit Hit.
No Improvement in 30 Generations

Best Fit Found at Generation 421
Fit Value = 4.230678e+01


Parameters at the Solution (value, gradient):

 X[ 1] :	5.656895e+00	G[ 1] :	1.774063e+01
 X[ 2] :	-6.534322e+00	G[ 2] :	-4.473903e-01
 X[ 3] :	6.877704e-01	G[ 3] :	-6.834089e-01
 X[ 4] :	5.467769e-01	G[ 4] :	-1.775774e+00
 X[ 5] :	-2.149052e+00	G[ 5] :	-2.170735e+01
 X[ 6] :	-3.893504e-01	G[ 6] :	-5.617223e+01
 X[ 7] :	7.203819e-02	G[ 7] :	1.146744e+02
 X[ 8] :	5.662580e-01	G[ 8] :	0.000000e+00
 X[ 9] :	1.740662e+00	G[ 9] :	8.071074e+01
 X[10] :	-7.724638e+00	G[10] :	-1.390066e+01
 X[11] :	9.045878e-01	G[11] :	1.344265e+00
 X[12] :	5.247284e-01	G[12] :	5.949332e+01
 X[13] :	-1.244724e+00	G[13] :	4.023601e+01
 X[14] :	-3.461685e-01	G[14] :	1.667676e+03
 X[15] :	3.870470e-01	G[15] :	-1.895165e+02
 X[16] :	1.941106e+00	G[16] :	4.004208e+01
 X[17] :	2.106121e+00	G[17] :	3.152157e+01
 X[18] :	9.727146e-01	G[18] :	0.000000e+00
 X[19] :	1.795138e+00	G[19] :	2.005074e+02
 X[20] :	-3.635496e+00	G[20] :	1.909344e+01
 X[21] :	7.922011e-01	G[21] :	3.054827e+01
 X[22] :	7.351241e-01	G[22] :	3.556713e+01
 X[23] :	-1.776021e+00	G[23] :	1.408034e+02
 X[24] :	-2.811830e-01	G[24] :	2.206455e+03
 X[25] :	3.441829e-01	G[25] :	8.791135e+02
 X[26] :	7.319743e-01	G[26] :	0.000000e+00
 X[27] :	-6.829757e+00	G[27] :	-1.127095e+02
 X[28] :	-3.023727e+00	G[28] :	-4.809519e+00
 X[29] :	2.334642e+00	G[29] :	2.701419e+00
 X[30] :	1.593285e-01	G[30] :	-1.023885e+01
 X[31] :	2.964882e+00	G[31] :	-8.462793e+01
 X[32] :	-3.625002e-01	G[32] :	-9.791735e+02
 X[33] :	1.720960e+00	G[33] :	5.164391e+01
 X[34] :	3.069959e+00	G[34] :	-1.616816e+01
 X[35] :	-5.168298e+00	G[35] :	1.065319e+00
 X[36] :	4.894117e-01	G[36] :	-5.701256e+00
 X[37] :	7.966906e-01	G[37] :	-7.082013e+00
 X[38] :	1.352425e-01	G[38] :	-1.811236e+01
 X[39] :	-4.014742e-01	G[39] :	4.203301e+01
 X[40] :	4.725082e-01	G[40] :	-1.021636e+02
 X[41] :	-3.259284e-01	G[41] :	0.000000e+00

Solution Found Generation 421
Number of Generations Run 433


Fri Mar  5 16:29:55 2004
Total run time : 1 hours 21 minutes and 16 seconds

(iter 1 ): LQD Results:
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.65689506  1.7406617  1.7951377 -6.8297569  3.0699594
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.53432165 -7.7246383 -3.6354955 -3.0237268 -5.1682985
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.68777038  0.9045878  0.7922011  2.3346415  0.4894117
SOE/SOE/SOE/SOE/SOE/NA                                                  0.54677689  0.5247284  0.7351241  0.1593285  0.7966906
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.14905230 -1.2447244 -1.7760211  2.9648822  0.1352425
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38935035 -0.3461685 -0.2811830 -0.3625002 -0.4014742
Age/Age/Age/Age/Age/NA                                                  0.07203819  0.3870470  0.3441829  1.7209604  0.4725082
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.56625801  1.9411059  0.7319743  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.1061212  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9727146  0.0000000  0.0000000  0.0000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD -0.3259284
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                         0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                  0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA               0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                    0.0000000
Age/Age/Age/Age/Age/NA                                                  0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000

(iter 1 ): LQD sigma: 42.30678 

WARNING: New fit is worse.  Stopping.
It may be a good idea to restart with a larger GENOUD population.
Warning message: 
NaNs produced in: sqrt(1 - Hdiag) 
> 
> summary(muljj, weights=TRUE);

Choice 1 : SLD Estimates and SE:
                 Est SE.Sand  SE.OPG SE.Hess t.val.Sand
(Intercept)    5.660  2.0500 0.32500 0.07220      2.770
New.Jobs      -6.520  2.2200 0.17700 0.06390     -2.940
Unemp          0.678  0.8080 0.13400 0.02890      0.839
SOE            0.544  0.3220 0.07220 0.01010      1.690
Church.Att    -2.160  0.3950 0.04980 0.01330     -5.470
Schooling     -0.389  0.0894 0.00753 0.00282     -4.360
Age            0.072  0.3460 0.05850 0.01240      0.208
Kwasn.SLD      0.571  0.0518 0.07030 0.00564     11.000
NA             0.000  0.0000 0.00000 0.00000        NaN
NA             0.000  0.0000 0.00000 0.00000        NaN


Choice 2 : PSL Estimates and SE:
                 Est SE.Sand SE.OPG SE.Hess t.val.Sand
(Intercept)    1.750  3.5600 0.3170 0.09010      0.492
New.Jobs      -7.710  3.2500 0.2810 0.07670     -2.370
Unemp          0.887  1.2200 0.1130 0.03430      0.725
SOE            0.520  0.3870 0.0753 0.01230      1.350
Church.Att    -1.230  0.5040 0.0589 0.01600     -2.450
Schooling     -0.346  0.1110 0.0116 0.00316     -3.130
Age            0.385  0.6180 0.0674 0.01710      0.622
Farmer         1.940  0.9310 0.0710 0.02130      2.080
Village        2.100  0.8120 0.0270 0.01630      2.590
Pawlak.PSL     0.975  0.0716 0.0703 0.00500     13.600


Choice 3 : UP Estimates and SE:
                 Est SE.Sand   SE.OPG SE.Hess t.val.Sand
(Intercept)    1.790  2.3300 4.34e-01 0.09640      0.767
New.Jobs      -3.630  1.4100 6.58e-01 0.08560     -2.580
Unemp          0.788  0.7970 2.44e-01 0.03900      0.990
SOE            0.734  0.2450 1.82e-01 0.01350      3.000
Church.Att    -1.780  0.4260 8.21e-02 0.01800     -4.170
Schooling     -0.281  0.0750 1.53e-02 0.00374     -3.750
Age            0.346  0.3590 8.30e-02 0.01650      0.965
Zioqkowska.UP  0.737  0.0792 2.24e+06 0.00535      9.310
NA             0.000  0.0000 0.00e+00 0.00000        NaN
NA             0.000  0.0000 0.00e+00 0.00000        NaN


Choice 4 : Catholic Estimates and SE:
                 Est SE.Sand SE.OPG SE.Hess t.val.Sand
(Intercept)   -6.880   5.080 0.2900 0.11200     -1.360
New.Jobs      -3.120   3.280 0.3710 0.09180     -0.950
Unemp          2.310   1.630 0.1280 0.04140      1.420
SOE            0.272   0.816 0.0482 0.01550      0.333
Church.Att     2.980   0.637 0.0816 0.01700      4.670
Schooling     -0.361   0.160 0.0122 0.00404     -2.250
Age            1.730   0.828 0.0557 0.01950      2.090
NA             0.000   0.000 0.0000 0.00000        NaN
NA             0.000   0.000 0.0000 0.00000        NaN
NA             0.000   0.000 0.0000 0.00000        NaN


Choice 5 : Other Estimates and SE:
                 Est SE.Sand  SE.OPG SE.Hess t.val.Sand
(Intercept)    3.060  2.0100 0.30100 0.06560      1.520
New.Jobs      -5.180  1.6800 0.29000 0.05640     -3.070
Unemp          0.502  0.9200 0.13900 0.02660      0.546
SOE            0.801  0.3010 0.09110 0.00913      2.660
Church.Att     0.138  0.4030 0.03930 0.01160      0.342
Schooling     -0.401  0.0864 0.00755 0.00260     -4.640
Age            0.474  0.3180 0.05640 0.01130      1.490
NA             0.000  0.0000 0.00000 0.00000        NaN
NA             0.000  0.0000 0.00000 0.00000        NaN
NA             0.000  0.0000 0.00000 0.00000        NaN


Choice 6 : UD.KLD Estimates and SE:
                 Est SE.Sand  SE.OPG SE.Hess t.val.Sand
Suchoka.UD    -0.331   0.237 3470000  0.0126       -1.4
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN
NA             0.000   0.000       0  0.0000        NaN


LQD sigma: 42.20445 
TANH sigma: 35.55004 

Number of Observations: 49 
Number of observations with at least one zero weight: 1 
Number of zero weights: 1 
TANH: weights
   name weights:SLD weights:PSL weights:UP weights:Catholic weights:Other
1     1    1.000000    1.000000          1         0.791814             1
2     2    1.000000    1.000000          1         1.000000             1
3     3    1.000000    0.834234          1         0.143349             1
4     4    1.000000    1.000000          1         1.000000             1
5     5    0.475479    1.000000          1         1.000000             1
6     6    1.000000    1.000000          1         1.000000             1
7     7    1.000000    1.000000          1         1.000000             1
8     8    1.000000    1.000000          1         1.000000             1
9     9    1.000000    1.000000          1         1.000000             1
10   10    1.000000    1.000000          1         1.000000             1
11   11    1.000000    1.000000          1         1.000000             1
12   12    1.000000    1.000000          1         1.000000             1
13   13    0.668669    1.000000          1         1.000000             1
14   14    1.000000    1.000000          1         1.000000             1
15   15    0.976474    0.976174          1         1.000000             1
16   16    1.000000    1.000000          1         1.000000             1
17   17    1.000000    1.000000          1         1.000000             1
18   18    1.000000    1.000000          1         1.000000             1
19   19    1.000000    0.846001          1         1.000000             1
20   20    1.000000    1.000000          1         1.000000             1
21   21    1.000000    1.000000          1         1.000000             1
22   22    1.000000    1.000000          1         1.000000             1
23   23    1.000000    1.000000          1         1.000000             1
24   24    1.000000    1.000000          1         0.559361             1
25   25    1.000000    0.847904          1         1.000000             1
26   26    1.000000    1.000000          1         1.000000             1
27   27    1.000000    1.000000          1         1.000000             1
28   28    1.000000    1.000000          1         1.000000             1
29   29    0.943129    1.000000          1         1.000000             1
30   30    1.000000    1.000000          1         1.000000             1
31   31    1.000000    1.000000          1         1.000000             1
32   32    1.000000    1.000000          1         1.000000             0
33   33    1.000000    1.000000          1         1.000000             1
34   34    1.000000    1.000000          1         1.000000             1
35   35    1.000000    1.000000          1         1.000000             1
36   36    1.000000    1.000000          1         1.000000             1
37   37    1.000000    1.000000          1         1.000000             1
38   38    1.000000    1.000000          1         1.000000             1
39   39    1.000000    1.000000          1         1.000000             1
40   40    1.000000    1.000000          1         1.000000             1
41   41    1.000000    1.000000          1         1.000000             1
42   42    1.000000    1.000000          1         1.000000             1
43   43    1.000000    1.000000          1         1.000000             1
44   44    1.000000    1.000000          1         1.000000             1
45   45    1.000000    1.000000          1         1.000000             1
46   46    1.000000    1.000000          1         1.000000             1
47   47    1.000000    1.000000          1         1.000000             1
48   48    1.000000    1.000000          1         1.000000             1
49   49    1.000000    1.000000          1         1.000000             1
> 
> save(muljj, file="jj.b.RData")
> proc.time()
[1] 4679.56  345.33 5039.33    0.00    0.00
> 
